Solitary Waves and Compactons 
in a class of Generalized Korteweg-DeVries 

Equations 

Fred Cooper 

Theoretical Division, Los Alamos National Laboratory, 
Los Alamos, JVM 87545 
and 

Physics Department, University of New Hampshire, 
Durham, NH 03824 

Harvey Shepard 

Physics Department, University of New Hampshire, 
Durham, NH 03824 

Pasquale Sodano 

Dipartimento di Fisica and Sezione, INFN Universita di Perugia, 
06100 Perugia, Italy 

February 9, 2008 



Abstract 

We study the class of generalized Korteweg-DeVries equations deriv- 
able from the Lagrangian: L(l,p) = J (h(p x <pt - j^pj + Oi(ip x ) p (ip xx ) 2 ') dx, 
where the usual fields u(x,t) of the generalized KdV equation are de- 
fined by u(x,t) = (p x (x,t). This class contains compactons, which 
are solitary waves with compact support, and when I = p + 2, these 
solutions have the feature that their width is independent of the am- 
plitude. We consider the Hamiltonian structure and integrability 
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properties of this class of KdV equations. We show that many of 
the properties of the solitary waves and compactons are easily ob- 
tained using a variational method based on the principle of least ac- 
tion. Using a class of trial variational functions of the form u(x, t) = 
A(t)ex.p — P(t) \x — q(t)\ 2n we find soliton-like solutions for all n, 
moving with fixed shape and constant velocity, c. We show that the 
velocity, mass, and energy of the variational travelling wave solutions 
are related by c = 2rEM~ 1 , where r = (p + I + 2)/(p + 6 — I), inde- 
pendent of n. 

PACS numbers: 03.40.Kf, 47.20.Ky, Nb, 52.35.Sb 



1 Introduction 

Recently, Rosenau and Hyman have shown that in a particular general- 
ization of the KdV equation, defined by parameters (m,n), namely 

K(n, m):u t + (u m ) x + (u n ) xxx = 0, (1) 

that a new form of solitary wave with compact support and width indepen- 
dent of amplitude exists. For their choice of generalized KdV equations the 
compactons with m = n < 3 had the form [cos(a^)] 2 ^ m ~ 1 \ where £ = x — ct 
and for m=2,3 they obtained: 

#(2,2) Mc = |cos 2 (e/4) 

#(3,3) u c = (!) 1/2 CO s(£/3). (2) 

Unlike the ordinary KdV equation, the generalized KdV equation consid- 
ered by Rosenau and Hyman was not derivable from a first order Lagrangian 
except for n — 1, and did not possess the usual conservation laws of en- 
ergy and mass that the KdV equation possessed. It is presumed that the 
generalized KdV equations found by the above authors are not completely 
integrable, but instead possess only a finite number of conservation laws. 
Because of this, we were led to consider a different generalization of the KdV 
equation based on a first order Lagrangian formulation. That is, we consider 
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L(l,p) = J (\p m - -^0^ + a( Vx y(ip xx ) 2 j dx. (3) 

This Lagrangian leads to a generalized sequence of KdV equations of the 
form: 

K*(l,p) : u t = u x u l ~ 2 + a (2u xxx ii p + ApvP~ x u x u xx + p{p - l)u p - 2 {u x f) (4) 
where 

u(x) = <f x {x). (5) 

These equations have the same terms as the equations considered by Rosenau 
and Hyman, but the relative weights of the terms are quite different leading 
to the possibility that the integrability properties might be different. [For 
the purposes of comparison it may be helpful to note that their set (m, n) 
corresponds to our (/ — 1, p + l).]The rest of the paper is organized follows: In 
section 2 we discuss some exact travelling wave solutions to (Q). In section 
3 we derive the conservation laws and discuss the Hamiltonian structure 
of these equations. In section 4 we apply the time dependent variational 
approach to obtaining approximate solitary wave solutions, and in section 5 
we compare the variational solutions to the exact ones. 

2 Exact solitary wave and compacton solu- 
tions 

If we assume a solution to (f|) in the form of a travelling wave: 

u{x,t)=f(Z) = f(x + ct), (6) 

one obtains for f: 

of = f'f- 2 + a (2f 'f + A P r l ff" + p(p - l)f p - 2 f' 3 ) . (7) 
Integrating twice we obtain: 
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We seek solutions where the integration constants, C\ and C2 are zero. This 
puts lower bounds on / and p: I > 1 and f"f p — » 0, / /2 / p_1 — > at edges 
where / — > 0. Then we obtain 

For finite /' at the edges, we must have p < 2,1 > p. 

Let us now look at some special examples. (Note that we have chosen 
signs so that all travelling waves have u > and move to the left.) The usual 
KdV equation has a — 1/2, I — 3, p — 0. For that case one has the well 
known soliton: 

u = (3c)sech 2 



3c/2£ 



(10) 



We define the "mass" M via 

/oo 
dx[u{x,t)] 2 . (11) 
-oo 

For this solution we find that we can express M and E in terms of c as 
follows: M = 24c 3 / 2 , E = f c 5 / 2 so that 

c = —EM' 1 = (M/24) 2/3 . (12) 
3 

The case / = p + 2 is the case relevant for compactons whose width is 
independent of the velocity c. For p = 1, a = 1/2 one obtains the compacton 
solution: 

Mi = 3ccos 2 (^/v^2), (13) 
where |f | < VZtt. One finds: M = f ttv^c 2 , E = f v^vrc 3 so that 

, ( 4M \ 1/2 

There is another compacton solution with p = 2, a = 3. 

u 2 = V6ccos(^/6) (15) 
with |£| < 37r. For this compacton, one finds M = 18irc, E = so that 

M 

c = AEM = . (16) 

187T 
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For the values, I = 3, p = 2 there is a compacton whose width depends 
on the velocity . Choosing a = 1/4 we find: 



u = 3c- (£ 2 )/6 



(17) 



on the interval 

If I < 3v^; (18) 

elsewhere it is zero. For this compacton one finds: M = ^v^c 5 / 2 , E = 
f v^c 7 / 2 so 



c= 11 jE ;M- 1 - , _. 
5 V144V/2/ 



2/5 



(19) 



Thus, apart from constants we find the same functional form for the 
compactons for our generalized KdV equations as those found by Rosenau 
and Hyman in their different generalization of the KdV equation. 

3 Conservation laws and canonical structure 

Equation (|J) can be written in canonical form displaying the same Poisson 
bracket structure as found for the KdV equation: 



5H 

u t = d x — = {u,H} 
ou 

where H is the Hamiltonian obtained from the Lagrangian (|3D, 



(20) 



H 



[(irip) — L] dx 



1(1-1) 

x 



u 



1(1-1] 



- a((f x ) p ((f xx ) 2 
dx. 



dx, 



- au p (u x ) 2 



(21) 
(22) 



By the usual arguments B this is consistent with a Poisson bracket struc- 



ture 



{u(x),u(y)} = d x 5(x - y). 



(23) 
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Let us now show that we have a system of equations which have exactly the 
same first three conservation laws as the ordinary KdV equation, namely the 
area, mass and energy. This is unlike the equations studied by Rosenau and 
Hyman that did not conserve the mass and energy, but instead had different 
conserved quantities. 
We have 

ut = dj¥- (24) 



5u 

so that the "area" under u(x, t) is conserved: 



lfu(x,t)dx = H 



Multiplying by u(x, t) we find: 



,u 



dt(-)=d x 



u 



+ a{(p-l)u p u 2 x + 2u p+1 u xx } 



which leads to the conservation of "mass" 

(1/2) J u 2 (x,t)dx = (1/2)M = H 1 



(25) 



(26) 



(27) 



For the KdV equation Hi was a second Hamiltonian under a second Poisson 
bracket structure. From Lagrange's equations we immediately get a third 
conservation law, the energy: 



H = 



u 



1(1-1) 



au p (u x Y 



dx = Ho. 



(28) 



The energy provided the first Poisson bracket structure: Considering the 
mass as a second Hamiltonian, the KdV equation has a second Poisson 
bracket structure using Hi. Assuming 



u t = {u,Hi} = J dy{u(x),u(y)}i 
one finds for the KdV equation that 



SHi 

5u(yY 



(29) 



{u(x),u(y)}i = (V + l -{Du + uD)j S(x - y) 



(30) 
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where D = d x . With this assumed Poisson bracket structure one again 
recovers the KdV equation. This Poisson bracket structure is identical to 
the Virasoro algebra with a specific central charge. This fact enables one 
to show that there is an infinite number of conservation laws in the KdV 
equation, and it is an exactly integrable system 

For the generalized KdV equations we find that we can write 

Ut = (a{D 2 u p D + Du p D 2 ) + \{Du l ~ 2 + u l ~ 2 D)\ u (31) 

so that there is a chance for a second Hamiltonian if the Jacobi identity is 
satisfied. One can postulate that the second Poisson bracket structure is 
given by 

{u(i),u(y)}i = (a(D 2 u p D + Du p D 2 ) + j (Du 1 ' 2 + u l - 2 D)^j 5{x - y). (32) 

So we need to show for what l,p this bracket structure obeys the Jacobi 
Identity, where the bracket is defined by : 

f°° SF 5G 

{F[u),G[u)} = ]_J-dy—{u(x)My)}i^y (33) 

One can show immediately that the Hamiltonians H\ and H2 commute using 
either Poisson Bracket structure ( p3|) or (|32|). 
We have that 

f°° 5H 2 SHi 
{H 2 [u],Hi[u}} = IJxdy—Mxluiy)}! — . (34) 

For the usual bracket structure (^3|) we can rewrite (0) as 



/oo 
dxu t (x 
-00 



Su(x) J 

-d t dxu 2 {x,t) = 0. (35) 
Z J —00 



2 

For the second bracket structure (|32D we have instead 

' 5H 2 \ 
8u{x)) 



/oo 
dxu t (x) 
-00 



7-oo Su(x) \5u(x) J 

- kL^imh - (36) 

Encouraged by this result we have attempted to repeat the induction 
proof of the existence of an infinite number of conservation laws, assuming 
as in the KdV equation that one has the conservation laws obey the recursion 
relations: 



(a(D 2 u p D + Du p D 2 ) + j(Du l ~ 2 + u l - 2 D)) = D^\. (37) 

V I J ou[x) ou[x) 

Starting with Hq defined by fl25|) we get the candidate Hamiltonian: 

f u l ~ l (x, t) , . , 

instead of (p27|). If we now ask if this is conserved by considering the equation 

^ = {H 1 ,H 2 } (39) 

using the first Poisson bracket structure, we find that the right hand side of 
is not a total divergence unless 1 = 3. For I = 3 one has : 

a(D 2 u p D + Du p D 2 ) + \{Du 1 - 2 + u l - 2 D)) = D-j^r (40) 

I ) bu[x) bu[x) 

However if we iterate one more time (with / = 3) we obtain: 

1,„ -A 5H. 

a 



{D 2 u p D + Du p D 2 ) + -(Du + uD)) = DF 3 (x) (41) 

3 / bu[x) 

and we find by explicit construction that F 3 (x) is not the variational deriva- 
tive of a local Hamiltonian unless p = 0. Thus this bi-Hamiltonian method 
of finding an infinite number of conservation laws only works for the original 
KdV equation. We surmise that (j3l|) is not a valid bracket structure and 
that 
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{H[u},{F[u},G[u}}} + {G[u},{H[u},F[u]}} + {F[u},{G[u},H[u}}} = (42) 



is not satisfied for the postulated second bracket. Thus we have not succeeded 
in showing that these new equations are exactly integrable, and we are in 
the same situation, in spite of having a first order Lagrangian, as for the 
generalized KdV equations of Rosenau and Hyman . 

We have not as yet performed numerical simulations of the scattering of 
our new compacton solutions. For Rosenau and Hyman such numerical ex- 
periments produced behavior very similar to but not exactly the same as that 
observed in completely integrable systems, namely, stability and preservation 
of shape. They find that elastic collisions are accompanied by the production 
of low amplitude compacton-anticompacton pairs ffl. 



4 Variational approach 

Our time-dependent variational approach for studying solitary waves is re- 
lated to Dirac's variational approach to the Schrodinger equation ||, |J. In 
our previous work || we introduced a post-Gaussian variational approx- 
imation, a continuous family of trial variational functions more general than 
Gaussians, which can still be treated analytically. Assuming a variational 
ansatz of the form u(x,t) = A(t) exp —/3(t) \x — q{t)\ 2ri ,we will extremize 
the effective action for the trial wave functional and determine the classical 
dynamics for the variational parameters. We will find that for all (l,p) the 
dynamics of the variational parameters lead to solitary waves moving with 
constant velocity and constant amplitude. For the special case of / = p + 2 we 
find immediately that the width of the soliton is independent of the amplitude 
and velocity. Correct functional relations between energy, mass, amplitude 
and velocity are obtained very quickly from the variational method, although 
one does not find that the I = p+2 variational solitons have compact support. 
We will find that most of the properties of the single "soliton" solutions to 
these equations can be obtained by using this very simple trial wave function 
ansatz and extremizing the action. 

The starting point for the variational calculation is the action 

Ldt, (43) 
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where L is given by (|3]). 

Just as we did in our study of the KdV equation we choose a trial wave 
function of the form: 



u v (x, t) = A(t) exp —j3(t)\x — q(t)\ 



\2n 



(44) 



where n is an arbitrary continuous, real parameter. 

The variational parameters have a simple interpretation in terms of ex- 
pectation values with respect to the "probability" 



P{x,t) 



[u v (x,t)]' 
M(t) 



where the mass M is defined as above 

M(t) = J [u v (x,t)] 2 dx. 

(Here we allow M to be a function of t, even though M is conserved) 
Since (x — q(t)) = 0, q(t) = (x). From (|46|) and ([P]) we have 



(45) 



(46) 



A(t) 



M l/2( 2 py/±r. 



2r (I + l) 



1/2- 



(47) 



The inverse width (3 is related to 

G 2n = (\x~q(t)\ 2n ) 



An(3 



(48) 



Following our approach in [ffl, we find that the action for the trial wave 
function ( fHP is given by: 



T(q, (3, M, n) = J (~Mq - C^P^^M 1 / 2 + C 2 (n)M 1+p / 2 /3 (p+4)/4n ) dt 
= jL!(q,q,M,P)dt, (49) 



where 
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CM 



C 2 (n) = 4«n(2)( p+2 )/ 4n (2 + p)^ 



v V4n 




(2-0/2 




Hi")'. 











-2 



r ( 2 "^) 



(50) 



[2r (X + l)]l +P /2- 

We eliminate the variable of constraint (3 (using 5T/5/3 = 0) and find 

(3 = [d{n)f n M 2n ^ +2 - l ^- p ~ & \ 



where 

d(n) 

From (|5TD we see that when 



(p + 4)C 2 (n) 
(Z-2)d(n) 

l=p + 2, 



(51) 
(52) 

(53) 



the width of the soliton (3 does not depend on M and thus is independent of 
the amplitude or velocity. This special case is precisely the case when the 
exact solution is a compacton. 

We now eliminate (3 in favor of M, and symmetrizing the Lagrangian (0) 
we obtain 1711: 



L 2 = - (qM-qM) — H(M), 



(54) 



where 



H(M) = (Cid (i - 2) - C 2 d {p+4) ) M r , (55) 
where r = (p + I + 2)/(p + 6 — I). Extremizing the action yields: 



M = 



M = const. 
(3 = const. 



(56) 



and 
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q = -2r (Cid (, ~ 2) - C 2 d^) M r ~\ (57) 
as well SIS db conserved energy 

E = (C X S 1 ^ - C 2 d (p+4) ) M r (58) 

Thus the velocity of the solitary wave is constant and can be related to the 
conserved energy via 



-c = —1rEM~ x . 



(59) 



This is precisely the form we obtained for the exact solutions. 

We have not yet extremized the action with respect to the variational 
parameter n which is equivalent to extremizing the energy with respect to n. 
We perform this extremization graphically for each value of Z,p. The explicit 
form of the trial wave function is: 



u 



v {x,t) = d[n,p,l]M 2 ^ p+6 -^2 1 / 4n 



2T 1 + 



exp 



2n 

_ d 4n M 2n(l-p-2)/( P +6-l), +ct _ 



-1/2 



Xq 



2)i 



(60) 



where d is given by (p2|). 

Now let us see how these trial wave functions compare with the exact 
answers for special cases. Since we explicitly know the M dependence of the 
answer, we can set M=l as our normalization for both the variational and 
exact solitons. 

First let us review the results for the KdV equation: Here the variational 
wave function is obtained by first setting p = 0,a = 1/2 and / = 3. One 
then extremizes the action in the constant parameter n. We find numerically 
that n=.877, which also extremizes the energy to be .035999M 5 / 3 and deter- 
mines the velocity to be .119995M 2 / 3 . In fi gure la and lb we compare this 
variational result to the exact soliton given by: 



(3c)sech 2 



3c/2(x + ct) 



(61) 



We see from fig lb that globally we achieve an accuracy of better than 1%. 
For this solution the velocity and energy are: 
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c = (M/24) 2/3 = .120187M 2/3 (62) 
3 

E = —Mc = .0360562M 5/3 . (63) 

10 V ; 

Thus we obtain the velocity (and the energy) accurate to 0.2% from the 

variational calculation. 

Next let us look at the compacton that is a segment of a parabola ([17|) . 

Here the variational wave function is obtained by setting p = 2, a = 1/4 

and I = 3. Mininimizing the action in the constant parameter n we find 

numerically that n=1.423, which also extremizes the energy to be .0803831 

and determines the velocity to be .225073 In figure 2a and 2b we compare 

this variational result to the exact soliton given by: 

(x + ctf 

u = 3c (64) 

6 

on the interval 

If I < 3^ (65) 

otherwise zero. We notice that the global accuracy is a few per cent except 
near the place where the true compacton goes to zero. For this compacton 
one finds: ^ 

c= (^_\ =.227006M 2 / 5 (66) 

and 

E = 5/UMc = .0810735M 7/5 (67) 

Thus we find that the velocity (and the energy) are determined to 0.8% 
accuracy. 

Next let us look at the compacton given by ([TBI). Here the variational 
wave function is obtained by setting p = l,a = 1/2 and I = 3. Extremizing 
the action in the constant parameter n we find numerically that n=1.154, 
which also extremizes the energy to be .054888 and determines the velocity 
to be .164666 In figure 3a and 3b we compare this variational result to the 
exact soliton given by: 

u x = 3ccos 2 (£/VV2), (68) 

where f < \^3n 
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We notice that the global accuracy is a few per cent except near the place 
where the true compacton goes to zero. For this compacton one finds: 



1 /2 

c= {^~7s) = - 165003Ml/2 ( 69 ) 

For the exact solitary wave one has the relationship: 

E = cM/2 = .0825017M (70) 

The approximate soliton had instead: 

E = cM/(2r) = cM/3 = .054888M (71) 

So for this compacton, the variational energy is wrong by a factor of 2/3 
although the velocity is correct to 0.2%. This is the only case where the 
variational method does not give the exact relationship between energy and 
velocity. However, we note that with a change in the sign of the second term 
in the expression for E, eq. (|58|), the variational energy becomes .0823332M, 
which is accurate to 0.2%, leaving the velocity and optimal n unchanged! 
(Such a sign change would follow from a factor (— l) p in C*2, which would 
leave all other results in this work unchanged.) 

Finally let us look at the compacton given by (15). Here the variational 
wave function is obtained by setting p = 2, a = 3 and I = 4. Extremizing the 
action in the constant parameter n we find numerically that n=1.283, which 
also extremizes the energy to be .00436284 and determines the velocity to be 
.017451 In figure 4a and 4b we compare this variational reult to the exact 
soliton given by: 

u 2 = Vfccos(£/6). (72) 

where £ < 3n 

We notice that the global accuracy is a few per cent except near the place 
where the true compacton goes to zero. For the exact compacton one finds: 

}l .0176839M (73) 



187T 

The variational estimates for c and E are accurate to 1.3%. For the exact 
and variational soliton one has the same relationship: 

M 2 

E = cM/A = — . (74) 
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FIGURE CAPTIONS 

Fig. la u v with n = .876 and u for M = 1 given by (61) as a function of x 
for the case: p — 0, a — 1/2, I — 3. 

Fig. lb u — u v as a function of x for p = 0, a = 1/2, I = 3. 

Fig. 2a u v with n = 1.423 and u for M — 1 given by (64) as a function of x 

for the case: p — 2, a — 1/4, I — 3. 

Fig. 2b u — u v as a function of x for p — 0, a — 1/2, / = 3. 

Fig. 3a u v with n = 1.155 and u for M = 1 as a function of x given by (68) 

for the case: p — 1, a — 1/2, I — 3. 

Fig. 3b u v — u as a function of x for p — 1, a — 1/2, I — 3. 

Fig. 4a u v with n = 1.283 and u for M — 1 given by (72) as a function of x 

for the case: p — 2, a — 3, Z = 4. 



15 



Fig. 4b u v — u as a function of x for p = 2, a = 3, / = 4. 
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